Splitting times of doubly quantized vortices in dilute Bose-Einstein condensates 
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Recently, the splitting of a topologically created doubly quantized vortex into two singly quan- 
tized vortices was experimentally investigated in dilute atomic cigar-shaped Bose-Einstein conden- 
sates [Y. Shin et al, Phys. Rev. Lett. 93, 160406 (2004)]. In particular, the dependency of the 
splitting time on the peak particle density was studied. We present results of theoretical simulations 
which closely mimic the experimental set-up. Contrary to previous theoretical studies, claiming that 
thermal excitations are the essential mechanism in initiating the splitting, we show that the combi- 
nation of gravitational sag and time dependency of the trapping potential alone suffices to split the 
doubly quantized vortex in time scales which are in good agreement with the experiments. We also 
study the dynamics of the resulting singly quantized vortices which typically intertwine — especially, 
a peculiar vortex chain structure appears for certain parameter values. 



Bose-Einstein condensation is associated with long- 
range phase coherence among the condensed particles. 
Due to this coherence, the dynamics of pure scalar con- 
densates can be accurately described using only a sin- 
gle complex valued order parameter. The fact that the 
condensate current density is proportional to the phase 
gradient of the order parameter implies vorticity to be 
quantized: vortices are topological defects in the order 
parameter field with a quantum number multiple k of 27r 
phase winding about the vortex line. 

Especially, properties of quantized vortices in gaseous 
atomic Bose-Einstein condensates (BECs) have been in- 
tensively investigated during the recent years, both ex- 
perimentally and theoretically, see Ref. m for a review. 
Since the angular momentum associated with a vortex 
is roughly proportional to the vorticity quantum num- 
ber K but the energy to its square, multiquantum vor- 
tices are typically energetically unfavourable. Conse- 
quently, the usual experimental methods to create vor- 
tices have yielded only singly quantized vortices or clus- 
ters of them However, using the topological phase 
engineering method originally suggested by Nakahara et 
al. 01, first two- and four-quantum vortices in dilute 
atomic BECs have been realized 0. This method uti- 
lizes the hyperfine spin degrees of freedom of the order 
parameter, but finally produces a scalar condensate con- 
taining a multiquantum vortex. 

Theoretical analysis has revealed that in addition to 
multiquantum vortices being energetically unfavourable 
in harmonic traps, they are also in general dynamically 
unstable ^ ^ . This dynamical instability implies that 
multiquantum vortices tend to split into singly quantized 
ones even in the absence of dissipation, i.e., in pure con- 
densates without noticeable thermal gas component. The 
splitting dynamics is an interesting issue, because the 
energy and angular momentum released from the multi- 
quantum vortex has to be redistributed in the system. 
The dynamics of the splitting of doubly quantized vor- 



tices created in the experiments reported in Ref. was 
theoretically studied in Ref. 0, and two major observa- 
tions were made: the initial dynamics in a cigar-shaped 
condensate can be modelled to some extent using only 
an effective two-dimensional local density analysis, and 
proper three-dimensional computations showed that the 
two vortices separating out from the doubly quantized 
vortex usually intertwine strongly as they split. 

Recently, the splitting of doubly quantized vortices was 
experimentally observed and the splitting time, i.e., the 
time interval between the creation of the vortex and the 
point when two separable vortex cores were observed, 
was measured as a function of the peak condensate den- 
sity 0- According to Ref. 0, the temperature was low 
enough for one to be able to neglect the effects of thermal 
atoms and dissipation. The experimental results verified 
that the doubly quantized vortex splits into two singly 
quantized ones. However, since only particle densities 
averaged over a short slice in the longitudal direction 
of the condensate were measured, no intertwining of the 
vortices was observed. 

In this Letter, we directly model the experiments re- 
ported in Ref. 01 , and compare the theoretical and exper- 
imental results. We solve the full three-dimensional dy- 
namics of the condensate using the Gross-Pitaevskii (GP) 
equation with time-dependent trapping potential com- 
bined with gravitational potential, closely corresponding 
to the experiments. From the condensate density pro- 
files, we determine the splitting time as a function of the 
peak condensate density, and analyse the effects of vor- 
tex intertwining. Contrary to previous theoretical results 
presented in Ref. we find that the gravitational and 
the time-dependent trapping potentials together break 
the rotational symmetry and initiate the splitting pro- 
cess strongly enough to alone yield splitting times in good 
agreement with experiments. Thus the effect of thermal 
excitations is not relevant in modelling the experiments. 

At low enough temperatures, the effect of the thermal 
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gas component can be neglected and the entire gaseous 
many-particle system of trapped atoms can be described 
by the condensate order parameter ^/;(r,t). Dynamics of 
this dilute condensate is governed by the GP equation 



2m 



V2 + l/(r,t)+g|V(r,t)| 



V-lr, t), 



(1) 

where V{r,t) is the external potential, m the mass of 
the atoms, and the strength of the interactions is gov- 
erned by the parameter g = Airh'^a/m expressed in 
terms of the s-wave scattering length a. The order pa- 
rameter is normalized according to the total number of 
atoms as / [■(/'pdr = TV. The stationary state solutions 
of the GP equation with eigenvalue are of the form 
■0(r, = e~'^'^*^^'^tp(r), and the excitation spectrum of 
these states can be solved from the Bogoliubov equa- 
tions. The excitation spectrum determines the stability 
properties of the state: Existence of modes with nega- 
tive energy but positive norm imply the corresponding 
stationary state to be energetically unstable, i.e., in the 
presence of dissipation and small perturbations the state 
will decay. Furthermore, existence of modes with non- 
real eigenfrequencies implies the state to be dynamically 
unstable, i.e., infinitesimal perturbations may grow ex- 
ponentially in time, and thus the stationary state may 
decay even in the absence of dissipation. 

Examples of these kind of states are the axisymmet- 
ric vortex states of the form ip{r) = e*^'^'^^-0(r, z), where 
(6,r,z) denote the cylindrical coordinates and k ^ 0. 
These states are energetically unstable in harmonic traps 
since vortices tend to spiral out of the condensate due to 
the buoyancy force. In addition, for quantum numbers 
|«;| > 2 there typically exists excitations with non-real 
eigenfrequencies which imply the corresponding states to 
be dynamically unstable and to split into singly quan- 
tized vortices. Especially, the dynamical instability of 
doubly quantized vortices has been demonstrated both 
theoretically 0, and experimentally 3 ■ 

In the experiments reported in Ref. 3|> Shin et al. 
created a doubly quantized vortex into a dilute BEG of 
^^Na atoms confined in a loffe-Pritchard trap by revers- 
ing the axial bias field linearly in 12 ms. The thermal 
gas fraction was not discernible in the experiments, and 
hence dissipation was presumably extremely weak. 

We have modelled these experiments by computing the 
time evolution of the condensate using the GP equa- 
tion. Starting from energy minimized axisymmetric dou- 
bly quantized vortex states corresponding to different val- 
ues of the density parameter auz — aj \ip{x,y,0)\'^dxdy 
(at the center of the condensate z — 0), we calculate 
the condensate dynamics in the time-dependent external 
potential corresponding to the loffe-Pritchard trap and 
gravity. Taking the gravity into account turns out to be 
crucial, since only the combined potential breaks the ro- 
tational symmetry about the vortex axis and is shown 



to initiate the vortex splitting process. In the compu- 
tations, we used finite-difference discretization. The ini- 
tial states were found using imaginary time integration, 
and the time evolution was calculated by a split operator 
method. 

In a loffe-Pritchard trap, the square of the total mag- 
netic field strength to the second order in the radial and 
axial coordinates is given by 



+ Bl{t), (2) 



where A determines the curvature of the initial axial field, 
C characterizes the strength of the field generated by the 
lofFe bars and B^ is the axial bias field. Including gravity, 
the external potential for the weak-field seeking state of 
spin-1 condensate reads 



V{T,t) = -gi^^iBBit) + Gmx, 



(3) 



where = ^1/4 and G is the gravitational acceleration. 
At t = 0, the term Bz{0) dominates in Eq. jSJ in the 
vicinity of the center of the trap, and hence the total 
potential is approximately harmonic in the condensate 
region. The initial potential is of the form 

^(r, 0) « imtj^z^ + imtj^[j/^ + (.T-a;o(0))^]+const. (4) 

According to the experiments, we choose the parameters 
A and C such that the axial and radial trapping frequen- 
cies are cu^ = 12 Hz and ujr = 220 Hz, respectively. In ad- 
dition, xo{t) = Gm\B,it)\/{gi^iiB[C'^ - Bz(t)A/2]) is the 
location of the minimum of (r, t) . The bias field B^ van- 
ishes at t = 6 ms, and hence the radial trapping potential 
is linear with respect to r. As Bz{t) passes through zero, 
the sign of A is changed in order to keep the axial field 
confining. After the reversal of at t = 12 ms, the 
potential assumes its initial form. 

The main results of our simulations are presented in 
Figs.n]-|Sl Figurc^shows the location xo{t) of the min- 
imum of y (r, t) and the computed location of the center 
of mass of the condensate as functions of time. Since 
the condensate slightly lags behind the time-dependent 
center of the potential, the rotational symmetry of the 
initial condensate state is broken. This perturbation also 
excites the dipole mode of the condensate, which is man- 
ifested by the oscillatory behavior of the center of mass 
after the potential has returned to its initial form. How- 
ever, the amplitude of the dipole oscillation is less than 
3% of the diameter of the condensate. 

In the effectively two-dimensional case, for which ojz = 
and the condensate is homogenous in the z-direction, 
the imaginary part of the frequency corresponding to the 
eigenmode responsible for the dynamical instability as- 
sumes a quasi-periodic form as function of auz 0, 0], 
and vanishes in some regions of the parameter aUz ■ Also 
in three-dimensional cigar-shaped condensates, regions 
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FIG. 1: Location xo{t) of the minimum of the total potential 
V{r,t) (dashed line) and the location Xcmit) of the center of 
mass of the condenstate (solid line). Due to inertia, the con- 
densate lags behind the center of the potential, which breaks 
the rotational symmetry of the system. This perturbation 
initiates the splitting of the doubly quantized vortex. 
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FIG. 2: (color online) Isosurfaces of condensate densities for 
(a) aUz = 2.56 at time t = 10.9 ms and (b) aUz — UA at 
t = 39.8 ms. The values defining the isosurfaces are of the 
order of 1% of the maximal density. 



where splitting is significantly faster or slower are clearly 
distinguishable. In Fig. E^a), the isosurface of a conden- 
sate with ariz — 2.56 is plotted at t = 10.9 ms. In this 
case, the splitting process is initiated along the whole 
length of the condensate, to be contrasted to the case 
aUz = 14.4 at t = 39.8 ms depicted in Fig.[2Ib), in which 
the splitting initiates only from the ends and center of 
the condensate. Since the precession frequency of the 
vortices depends on their distance, intertwining is ob- 
served in cases which contain such stable regions. Note 
also the distinctive surface mode excitations due to the 
additional energy and angular momentum released in the 
splitting process. 

We determine the splitting times T from the den- 
sity profiles n{x,y) = /j^ig^^^^ \^{x,y, z)\'^dz, mimicking 
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FIG. 3: Density profiles n{x,y) = ^_^^l^'^^\il){x,y,z)\'^dz of 
a condensate with an^ = 2.56 at three different times t. 
The splitting time of the doubly quantized vortex is deter- 
mined by comparing the separation of the local minima in 
the density profiles to the diameter of a single quantum vor- 
tex represented by the white circle corresponding to a diam- 
eter of 1.8 /im in (a). The field of view in these figures is 
13.5 /im X 13.5 /im. 

the tomographic imaging technique used in the experi- 
ments 0- Roughly, when two minima are observed in 
the density profile with their separation exceeding the 
diameter of a singly quantized vortex, the doubly quan- 
tized vortex was considered split in the experiments. We 
define the diameter of a single quantum vortex as twice 
the distance from the vortex center to 75% of the maxi- 
mum value in the density profile n{x,y) of the stationary 
single quantum vortex state. In Fig. |31 the density pro- 
file fi{x,y) is plotted at t = 10.9,14.5 and 18.1 ms for 
aUz = 2.56. From the isosurface plot in Fig. |5fa), one 
might argue that the vortex has already split at t — 10.9 
ms, but from the density plot in Fig. 01 a) no such con- 
clusions can be made. In Fig. O^b), the distance of the 
local density minima already exceeds the single quantum 
vortex diameter and a detailed analysis yields T = 13.7 
ms for the splitting time. 

Figure 0] shows the splitting times for values of aUz 
ranging from to 20 corresponding up to = 8 x 10^ 
^^Na atoms. A finite splitting time T is obtained also 
in the noninteracting case, which concurs with the ob- 
servation that persistent stationary currents exist in the 
noninteracting case only if the spectrum of the Hamilto- 
nian is degenerate For aUz < 3, T decreases with 

increasing interaction strength. This effect is also visible 
in the recent theoretical study of Ref. j^, but it is not 
clearly observable in the experimental results 0. The 
splitting time attains its minimum at aUz ~ 3, which is 
in good agreement with the Bogoliubov eigenvalue spec- 
trum analysis accomphshed in Refs. ^'^^ '^'^z > 10, 
T saturates to a finite value, which is in fair agreement 
with the experimental results although the satura- 
tion value is roughly 20% larger in the experiments. The 
overall form of the computed splitting times as a func- 
tion of aUz agrees well with the experimental data, and 
the agreement is surprisingly good even for the quantita- 
tive results. Note that the computed splitting times are 
generally somewhat smaller than the measured ones, im- 
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FIG. 4: The splitting time T of a doubly quantized vortex 
as a function of ariz. The splitting time decreases rapidly as 
interactions are introduced in the system attaining its mini- 
mum value around auz « 3. As the interaction strength is 
increased further, T increases due to stable region near the 
center of the condensate. The splitting time saturates roughly 
to T ~ 55 ms for interaction strengths ariz > 10. 
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before imaging, and in our results possible changes in the 
relative vortex separation during this expansion period 
are neglected. 

For values near an^ ~ 10, where there is a stable re- 
gion at the center of the condensate, a peculiar linked 
chain-type structure for the intertwined vortices was ob- 
served, see Fig.lS^a). The ^-integrated density profile ob- 
served with the tomographic imaging technique is shown 
in Fig.E^b). In this case, there are four distinct minima 
in the density profile instead of the usual two. This effect 
could explain the observation interpreted as crossing of 
the vortex lines in the experiments. 

In conclusion, we have calculated the splitting times of 
doubly quantized vortices as a function of peak particle 
density by solving the time-dependent Gross-Pitaevskii 
equation numerically, mimicking closely the experiments 
reported in Ref. 0. The results are in good agreement 
with the experimental data, confirming that thermal fluc- 
tuations are not required to explain the measured split- 
ting times. Instead, the main impetus for the rotational 
symmetry breaking and splitting is the combination of 
the gravitational and time-dependent trapping poten- 
tials. 

The authors thank CSC, the Finnish IT Center for Sci- 
ence, for computational resources. J. H. and M. M. ac- 
knowledge the Finnish Cultural Foundation for financial 
support. Y. Shin is appreciated for helpful discussions. 



FIG. 5: (color online) (a) Isosurface plot of the vortex chain 
structure formed in a condensate with ariz = 11.0 at t = 
54.3 ms. The value defining the isosurface is of the order of 
10% of the maximal density, (b) The corresponding density 
profile n{x,y) shows four minima instead of two. The field of 
view in (b) is 6.8 /xm x 6.8 pm. 



plying that the time dependency of the external potential 
alone gives sufficiently strong impetus for the splitting, 
and the additional contribution of thermal fluctuations 
is not needed. 

A few points should be taken into account when com- 
paring the computed splitting times to experimental val- 
ues. Firstly, our simulations neglect the multicomponent 
nature of the condensate when the multiquantum vortex 
is created by reversing the bias field B^. In our sim- 
ulations, we begin with the multiquantum vortex state 
already at t = 0. Effectively, this difference can roughly 
be assumed to shorten the computed lifetimes of the mul- 
tiquantum vortices by half of the reversing time, i.e., by 
6 ms. However, this difference is presumably more than 
compensated by the fact that in the simulations we start 
counting the splitting time from the beginning of the per- 
turbation, whereas in the experiments the time is counted 
from the end of the perturbation. Secondly, in the exper- 
iments the condensate is let to expand freely for 15 ms 
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